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Seismic,  depth,  depth  phase,  echo  detection 


#<>•  AO?TKaCT  (Ceniinui  on  rt'Vtif  if  n*  < ;J  ir'rntifv  by  block  nnnl’er) 

^Several  modi f icati ons  were  made  in  the  seismic  source  depth  deter 
mination  program.  Ihese  m.od  i 1 icat  i ons  were  uirectcil  towards 
improving  tlie  proprr.ti ' s use  of  depth  pliase  information,  and 
increasing  tlie  i nterpreta!' i 1 i t y of  tlie  final  depth  plot  outinit. 
Two  new  events  were  run  throu;’,ii  tiiis  modified  version  of  tlie 
program,  with  successful  deptii  determinations  being  made  in  hotli 
case  s . 
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INTRODUCTION  AND  SUMMARY 


In  the  second  quarter  of  this  contract,  a final  set  of  mod- 
ifications was  made  to  tlie  seismic  source  depth  determination 
program,  and  the  resulting  production  version  was  run  on  five  new 
events.  The  production  version  was  decided  on  after  development 
and  testing  of  a new  data  windowing  teclinique  designed  to  reducc 
the  occurrence  of  unwanted  cepstrum  peaks.  Although  this  new 
technique  produced  the  desired  results  in  some  cases,  generally 
it  did  not  work  as  well  as  tlic  older  windowing  method;  and,  con- 
sequently, the  older  method  was  the  one  used  in  the  produ''*ion 
program.  Changes  incorporated  in  tlie  production  version  include 
the  capability  of  obtaining  depth  plots  using  several  different 
cepstrum  computation  window’  lengths,  and  a peak  significance  level 
calculation  designed  to  assist  in  tlie  interpretation  of  tlie  dept  i 
plots.  Results  obtained  from  the  five  new  events  are: 


Event 

Date 

No.  of  Stations 
Used 

Depth  (km) 

3/27/75 

4 

17.5+1.5 

6/30/75 

5 

12.. 5+2.0 

3/7/75 

2 

19.5  + 1 . 5 

4/28/75 

2 

not  determined 

6/1/75 

3 

13.5+2.0 

MAJOR  ACCOMPLISliMHNTS 


PROGRAM  DIATI.OPMHN'T 


Nc's  Drita  V,’i  lidov;  i n»  Tcchni 

Two  new  techniques  for  positioning  the  seisinograiu  windows 
used  in  cepstruin  computations  were  developed  and  tested.  The 
principal  goals  of  tliese  new  methods  were: 

1 . improvement  of  cepstrum  peaks  due  to  depth  phase 
correlations;  and 

2.  reduction  of  spurious  cepstrum  peaks  due  to 
correlations  between  different  primary  phases. 

The  first  of  tl^ese  new  tecJiniques  consists  simply  of  allowing 
windo\;s  to  ovealap  instead  of  using  only  adjacent,  non -overlapping 
seismogram  windows.  This  should  make  it  more  likely  that  a window 
will  exist  that  contains  both  a primary  arrival  and  its  correspond- 
ing depth  phases,  resulting  in  a cepstrum  tliat  contains  tlie  desired 
complete  set  of  depth  phase  peaks.  The  second  method  involves 
defining  different  cepstrum  computation  windows  for  each  primary 
phase,  station,  and  trial  depth.  Each  windov>'  starts  at  the  expected 
primary  phase  arrival  time  and  ends  after  tl;e  predicted  arrival  of 
the  corresponding  s-phase.  This  metliod  guarantees  that  depth  phases 
will  always  fall  in  the  same  window  as  their  corresponding  primary 
phase,  and  should  reduce  the  occurrence  of  spurious  cepstrum  peaks 
by  always  using  the  shortest  window  necessary  to  include  tlie  com- 
plete set  of  expected  depth  phases. 


r 


} 

I 


The  now  windowing  techniques  were  compared  to  the  original 
windowing  method  using  five  different  events.  TIic  results  indi- 
cated tliut  t’ne  constant  length,  non-overlapping  windows  that  were 
originally  used  worked  best.  Although  tJic  arrival  - time-control  1 eJ 
windows  gave  good  results  in  some  cases,  and  almost  always  showed 
some  evidence  of  the  correct  dcptii  peak,  the  results  from  the 
constant  windows  were  more  generally  dependable.  This  somewhat 
unexpected  o'.’.tcome  may  be  due  to  scaling  inconsistencies  result- 
ing from  the  use  of  variable  cepstrum  computation  window  lengths. 
The  basic  idea  bcliind  tliis  windowing  method,  liowever,  seems  to  be 
a useful  one,  and  it  may  merit  further  investigation  at  some 
future  time. 

Implementation  of  Production  Program 


Based  on  the  results  of  the  investigation  of  seismogram  window- 
ing techniques,  tlie  original  version  of  the  depth  determination 
program  was  modified  to  serve  as  the  production  version.  Two  ncv.- 
capabilities  have  l^cen  added  in  the  production  program. 

1.  Processing  of  seismograms  at  several  different 
cepstrum  computation  window  lengths. 

2.  Determination  of  a significance  level  to  be  used 
in  interpreting  final  depth  plot  output. 

Previous  experience  has  indicated  that  these  two  program 
features  should  be  particularly  useful  in  production  processing  of 
large  numbers  of  events.  Depth  plot  output  has  been  found  to  be 
very  sensitive  to  cepstrum  computation  window  length,  and  the  opti- 
mum window  length  varies  unpredic tal)l y between  events.  Also,  it 
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is  often  not  clear  which,  if  any,  of  the  peaks  appearing  on  the 
depth  plot  output  are  meaningful.  A quantitative  way  of  estimat- 
ing the  significance  of  these  peaks  would  be  very  useful  to  tiie 
analyst. 

A teclinique  lias  been  implemented  for  estimating  tlie  signifi- 
cance of  depth  plot  amplitudes  and  is  based  on  the  fact  that  these 
amplitudes  relate  to  the  degree  cepstrum  amplitudes  sum  when 
selected  at  lag  times  expected  for  a given  source  depth  and  epi - 
center  distance.  Therefore,  one  way  to  obtain  a measure  of  depth 
plot  amplitudes  which  would  result  from  the  analysis  of  seismic 
data  containing  the  same  spectral  characteristics  but  not  con- 
taining depth  phase  information,  is  to  randomize  the  travel  time 
information  and  reprocess  a given  event.  In  this  way,  one  can 
simulate  a depth  plot  resulting  from  a seismic  recording  in  which 
the  cepstrum  peaks  did  not  result  from  depth  phase  information 
and  do  not  constructively  add  for  travel  times  corresponding  to 
any  given  source  do]ith.  For  example,  when  ana.lyzing  data  recorded 
at  A=30°  for  a trial  depth  of  20  km,  the  program  would  normally 
add  cepstrum  anplitudes  at  delays  of  S.6  seconds  (pP-P) , 8.3 
seconds  (sP-P),  and  2.7  seconds  (sP-pP).  In  computing  the  signi- 
ficance level,  a depth  plot  amplitude  at  20  km  will  be  obtained 
by  summing  cepstrum  values  at  three  random  times  sucli  as  8.1,  5.4 
and  4.7  seconds.  A random  pick  deptli  plot,  like  the  one  shown 
in  Figure  la,  is  constructed  Iiy  making  random  picks  like  tliis  for 
each  dcptli  being  considered.  A histogram  of  the  random  depth  plot 
amplitudes  is  then  constructed  (Figure  lb),  and  a cumulative  dis- 
tribution function  (Figure  Ic)  is  calculated  from  the  liistogram. 
Finally,  the  amplitudes  corresponding  to  tlie  80  percent  and  95 
percent  levels  on  the  cumulative  distribution  are  determined  and 
plotted  on  the  true  depth  plot  (Figure  Id).  V.’hat  is  indicated  by 


tliesc  levels  is  the  prolKil)  i 1 i ty  th;it  the  resulting  peaks  could 
not  have  been  produced  by  cliance  alignment  of  cepstrum  peaks. 

Tlius , pealcs  above  the  95  percent  level  have  a 95  percent  confid- 
ence that  tliey  did  not  occur  by  cliance  cepstrum  peak  alignment. 

However,  to  establish  the  percentage  confidence  of  a peak  above 

the  95  percent  level  relating  to  the  correct  event  dcptli,  results  j 

from  a large  set  of  events  must  be  compiled.  ^ 

t 

Nf.w  Evr.irfs 

I-ive  new  events  were  run  during  this  quarter,  witli  deptli  j 

estimates  being  successfully  obtained  for  four  of  them. 

3/27/75  Event  n'urkcyi 

Seismograms  from  the  CA’^ent  dated  3/’27/75  (Turkey)  are  shown 
in  Figure  2.  The  best  results  were  obtained  using  the  first  6.4  ^ 
seconds  of  data  from  all  four  stations,  with  a 12.8  second  cep- 
strum  computation  window;  the  resulting  composite  dcptli  ]ilot  is 
shown  in  Figure  3.  Interpretation  was  quite  easy  for  this  event, 
with  a single  dominant  peak,  well  aboA'o  the  95  percent  significance 

■ 

level,  appearing  at  17.5  km.  These  excellent  results  are  due  to  j 

the  strength  of  the  depth  phase  arrivals  for  this  CA’^ent  - pP  can  | 

be  visually  detected  on  three  of  the  four  seismograms  in  Figure  2.  | 

! 

1 

! 

6/30/75  liven t (Montana)  j 

Seismograms  from  the  oA’^ent  dated  6/30/7  5 (Montana)  are  shown  j 

in  Figure  4.  The  best  results  were  olitained  using  the  first  51.2  j 

seconds  of  data  from  all  five  stations,  with  a 12.8  second  cepstrum  I 

computation  window;  the  resulting  composite  depth  plot  is  shown  | 


-6- 


Turkey  Event  {Z/21flS)  Scismograr.s 
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Montana  Event  (6/50/75)  Seismograms 


in  ]-igu!'c  5.  This  event  is  not  as  clear  cut  as  tlie  previous  one, 
hut  tiie  peak  at  12.S  km  extc'iids  far  enough  above  the  ‘>S  percent 
significance  levc'l  to  be  picked  as  tlie  correct  dc])tli.  A likely 
explanation  for  tlie  poorer  finality  of  the  depth  jilots  ol)tained 
for  tliis  event  can  be  found  by  looking  at  the  seismograms.  In 
contrast  uitli  t!ie  previous  event,  no  dc'pth  phases  can  i)C  seen  at 
any  of  tlie  five  stations. 

>/! /IS  I'.vent  C Irani 

Til  i s event  jiresented  the  most  difficult  interpretation  pro- 
blem of  all  the  events  in  this  set,  and  i\ill  be  discussed  in  more 
detail  than  the  others.  Seismograms  for  this  event  arc  shown  in 
Figure  6,  ami  Figures  7-P  are  the  i’.  Pci’,  and  composite  dejitli 
plots  ol'tained  using  .SI.  2 seconds  of  data  from  both  stations,  with 
a 12.8  second  cepstrum  computation  window.  From  the  composite 
depth  plot  alone,  it  aiipears  that  tlic  pc-ak  at  8 km  should  be 
interpreted  as  indicating  the  correct  depth.  However,  tlierc  are  ■ 
some  prclilems  witli  this  int  erjirctat  i on . First,  the  liciglit  of  tlie 
8 km  peak  above  the  95  percent  significance  level  probably  is  not 
very  meaningful,  since,  with  only  two  stations  being  used,  the 
statistical  basis  of  the  significance  level  calculation  is  not 
very  good.  Second,  the  P-phasc  depth  plot  shows  no  single  domiimnt 
peak,  and  the  PeP  plot  shows  that  the  8 km  peak  on  the  composite 
plot  comes  from  the  reinforcement  of  the  strong  8 km  PeP  peak  and 
tlic  much  weaker  8 km  peak  seen  on  P.  This  reinforcement  is  not, 
in  this  case,  indicative  of  a depth-consistent  cepstrum  peak.  At 
the  large  A's  being  considered,  P and  PeP  arc  being  picked  from 
the  same  cejistrums  and  have  the  same  delay  times,  so  any  peaks 
occurring  in  the  first  cepstrum  comjuitalion  window  (used  to  generate 
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No.  of  stations  = S 

Ccpstruni  windov:n  length  = 12.8  see 


Montana  livent  (0/80/78],  (iomposite  I'eplli  I’lol 


Iran  Event  (3/7/75)  Scismograns 


No.  of  stations 


Iran  Hvcnt  f5/'^/73')  , P Pcplli  Plot 


I’cF  l)c'i)t!i  i’lut 


.\o.  ot  stations 
Copstruin  uindoK  length 
Total  Icngtli 


lompos 1 to 


both  P ontl  PcP  dcptli  plots]  wilt  be  enl):inccd  over  pcoks  from  later 
cepstrums  (used  only  for  P plot).  Conseciuent  1 y , the  composite 
will  contain  peaks  tliat  liave  l)een  unfairly  reinforced,  and  only 
the  P de])tli  plot  v.’ill  be  useful  for  interpretation. 

A better  result  was  obtained  \;hen  tlie  same  data  was  processed 
Avith  ;i  sc'cond  cepstrum  computation  window.  The  P depth  plot, 

shown  in  figure  10,  lias  a single  dominant  peak  at  19.5  km,  which 
can  justifiably  be  interpreted  as  indicating  the  correct  dejitli. 

Tlic  success  of  the  depth  determination  ])rogram  for  this  event  is 
cncoui'aging  in  view  of  the  small  numlier  of  stations  used  and  tlie 
lack  of  visually  detectable  dejith  phases  on  the  seismograms. 


4/28/75  hvent  (Kasmir) 


This  event  is  tlie  only  one  of  the  five  events  processed  for 
which  no  depth  determination  could  be  made.  The  seismograms  are 
shown  in  Figure  11.  Like  the  Iran  evenv.,  both  stations  have 
large  A's,  so  only  the  P depth  plots  are  likely  to  be  useful. 

The  P depth  plot  obtained  using  51.2  seconds  of  data  from  both 
stations,  with  a cepstrum  computation  window  of  12.8  seconds  is 
shown  in  Figure  12.  A peak  is  present  at  17  km,  but  its  lov;  signi- 
ficance level  makes  it  unacceptable.  Figure  1.5  is  the  P plot 
from  the  same  data,  with  a 25.6  second  cepstrum  computation  window. 
There  is  a strong  peak  present  at  4 km,  but  this  depth  represents 
a pP-P  time  of  about  one  second,  which  is  too  short  to  lie  detected 
using  present  methods. 
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Kasmir  Event  (4/28/75),  P Depth  Plot 


6/1/ 7 5 Eren  t_  [ Mo  i a v c l ie  sort) 

Seisiuograms  for  tliis  event  arc  sliown  in  Figure  14.  Best 
results  v;ere  obtained  using  a 12.8  second  cepstruiii  computation 
windoh'  on  76. S seconds  of  data  froin  all  tliree  stations;  the  coin- 
po.site  depth  plot  is  shorn  in  I’igure  15.  The  doninanr  peak  at 
13.5  km  is  well  aboA'e  the  95  percent  sign  i f ic.ance  level,  and  ras 
interpreted  as  indicating  the  correcl  de])th. 

Coiicl  us  ions  From  Analys  is  of  \cv.^  F.ven  ts 

Several  conclusions  can  be  drawn  from  the  depth  determination 
program  results  for  these  five  events. 

• Results  from  tlie  Turkey  and  MiJ.ire  Desert  events 
sliow  tliat  source  deptlis  from  tlie  depth  determina- 
tion program  agrees  wi  tli  tlios('  obtained  by  observ- 
ing pP  on  the  seismograms  wlien  pP  is  visible. 

• Results  from  the  Montana  and  Tran  events 
demonstrate  that  source  depth  estimates  can 

be  obtained  even  wlicn  no  depth  phases  arc  visible. 

• Evidence  to  date  indicates  that  analysis  involv- 
ing fewer  tlian  three  stations  significantly 
reduces  the  proinability  of  obtaining  a depth 
estimate . 

• The  computation  of  significance  levels  appears 
to  be  a very  useful  technicpie  and  should  be 
investigated  further. 
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Firm  RE  PLANS 


In  the  remainder  of  this  contract,  we  plan  to  process  tiie 
15  station  Andreanof  Island  event  (11/22/65)  in  several  subsets 
to  furtlior  investigate  the  effects  of  the  number  of  statio;is  used. 
Any  new  events  tliat  are  provided  will  also  he  run.  Results  from 
all  tlie  events  tliat  liave  been  processed  will  be  analyzed  for  tlie 
effects  of  number  of  stations,  seismogram  signal - to-noise  ratio, 
and  station  I's. 


l7i  the  area  of  program  development,  present  results  indicate 
the  need  for  future  v,-ork  in  two  areas.  First,  improvements  in 
the  determination  of  peak  significance  may  be  possible.  For 
e.xample,  meaningful  estimates  of  significance  levels  approaching 
100  percent  m.a}  be  possible  if  some  kind  of  curve  fitting  is  done 
on  the  random  pick  cumulative  distrilnit; on  function.  Second, 
techniques  for  eliminating  cepstrum  peaks  due  to  non-deptli  phase 
spectrum  structure  need  to  be  e.x'plorcd.  This  may  eliminate  the 
shallow'  depth  "noise"  observed  on  many  depth  plots,  and  aid  in 
the  detection  cf  shallow  depth  eartliquakes . 
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